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We present a full characterization of the phase transition in U(l) scalar field theory and of the as- 
sociated vortex string thermodynamics in 3D. We show that phase transitions in the string densities 
exist and measure their critical exponents, both for the long string and the short loops. Evidence 
for a natural separation between these two string populations is presented. In particular our results 
strongly indicate that an infinite string population will only exist above the critical temperature. 
Canonical initial conditions for cosmic string evolution are show to correspond to the infinite tem- 
perature limit of the theory. 
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Topological defects appear in a great variety of sys- 
tems from condensed matter laboratory experiments to 
the early Universe. Their importance in phase transi- 
tions in the laboratory is know to be fundamental and 
their presence in the early Universe may be the key to 
many of the unsolved questions in standard cosmology. 

However, in spite of the universal relevance of topolog- 
ical defects, much about the fundamental description of 
their formation and evolution remains qualitative. This 
is a reflection of the complexities involved in first prin- 
ciple studies, owing to their nature as non-perturbative 
excitations of quantum field theories. 

Particularly interesting for cosmology are string-like 
topological defects [Q . Motivated by the study of their 
creation in the early Universe [Q] , a variety of experiments 
has been recently developed with the aim of studying 
vortex string formation in liquid crystals Sj, superfluid 
4 He H and 3 He || systems. Their new results permit 
us to test with unprecedented precision theoretical ideas 
about defect formation and evolution. 

From a theoretical standpoint, cosmic strings and 
other topological defects, have been traditionally thought 
to be produced at phase transitions in the early Universe 
as a result of the formation of correlated domains 
More recently a refinement of this scenario has been 
gaining support, that most topological defects existing 
below the phase transition are, in fact, survivors of a 
population of unstable defects existing above the phase 
transition as non-linear excitations of the fields. In order 
to exist at low energies these field configurations need 
to be frozen in by some non-equilibrium process, which 
is context dependent. In a continuous phase transition, 
like the superfluid transition in 4 He, this process may in 
addition be helped by the critical slowing down in the 
response of the fields close to the critical point. An es- 
timate of the defect density hence produced is usually 
achieved given the correlation length of the field at the 
time of freeze out. By assigning random phases to dif- 
ferent domains and connecting by the path of minimum 
phase gradient, a network of defects can in principle be 



constructed. This picture generates order of magnitude 
estimates that successfully predicted the defect densities 
observed in the Helium experiments. 

In practice higher correlations in the phase structure of 
the fields exist and can change the picture considerably. 
As an illustration of the shortcomings of simply using 
the correlation length to predict produced defect densi- 
ties note that similar densities would be estimated if the 
defect network were frozen in above or below the critical 
point, provided that the correlation length were chosen 
to be the same in both cases. This is clearly not what 
is observed experimentally jj| , suggesting that the string 
network must be very different in two circumstances were 
the domain structure can be expected to coincide. 

In order to match theory to future, more accurate mea- 
surements, as well as for the sake of theoretical under- 
standing it is highly desirable to generate more detailed 
predictions for the string density formed at phase tran- 
sitions. In cosmology and in the laboratory, equilibrium 
is the natural starting point. If we know the statistics of 
strings at any given temperature the reference to the do- 
main structure in the fields becomes obsolete, while the 
role of non-equilibrium in freezing in the string network, 
at least in some scales, can be expected to hold as usual. 

In this Letter we present the full characterization of the 
behavior of the string densities with temperature in the 
U(l) scalar field theory. In doing so we go beyond exist- 
ing studies in the XY (see eg. g]) model, which belongs 
to the same universality class, especially in determining 
the behavior of infinite strings, measuring critical expo- 
nents associated with string densities and confronting our 
findings with cosmological scenarios of defect formation. 

We first determine the field thermodynamics and lo- 
cate the critical point. In order to drive the system to 
equilibrium we evolved the fields stochastically according 
to 
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where {i,j} <E 1,2. The stochastic variable £i(x, t) is 
taken to be Gaussian, characterized by 

&(x,t)) = 0, (2) 

&(x,t)£,-(x',0> = ^(x - x')«5(i - *')%• (3) 

These relations ensure that detailed balance applies in 
equilibrium, which will result for large times at tem- 
perature We evolved the system Eq. (0-||) using 
a staggered-leapfrog method, on a lattice of size 100 3 . 
We chose the lattice spacing to be SI = 0.5. We verified 
that the correlation length was always resolved by at least 
4 points in each linear dimension, apart from the cases 
when /3 < 1. Around the critical point all physical scales 
become much larger than the lattice spacing. Universal 
quantities can then be measured and shown explicitly to 
be independent of this ultraviolet cutoff. Away from the 
critical point, at high or low temperatures, universality 
is lost in general. There, the results loose some of their 
physical meaning but retain model properties, that are 
interesting from a theoretical point of view. This will be 
the case of the infinite temperature limit of the theory. 

An equivalent equilibrium state could have been ob- 
tained using lattice Monte Carlo methods. Also, from the 
equilibrium point of view, the second order time deriva- 
tive is redundant. However, this work is intended to be 
the precursor to non-equilibrium studies of the relativis- 
tic field theory, where is is more natural to use a Langevin 
method such as ours, and the field equations are natu- 
rally second order in time derivatives. 




FIG. 1. Phase probability distributions around the critical 
temperature for (3 between 1.1 and 2.9 (steps of 0.1). For 
low temperatures the distribution becomes peaked around a 
randomly chosen phase. As the temperature increases the 
distribution becomes flat and the U(l) symmetry is restored. 

We decide when the system has reached equilibrium 
by monitoring several quantities associated with differ- 
ent length scales. We measure the kinetic energy of the 
system and check for equipartition. This characterizes 
the smallest scales in the sample. We monitor the string 
densities, associated with intermediate length scales and 
we measure the phase distribution functions across the 



whole volume. By making sure that these three quanti- 
ties have reached asymptotic average values we conclude 
that the system is in equilibrium. The temperature of the 
bath coincides with that measured using equipartition. 

The transition between a system displaying a prefered 
random direction in its manifold of field minima at low 
temperature and a flat phase distribution, characteristic 
of the U(l) symmetry is shown in Fig. [I]. 

Although appealing the phase probability plots of 
Fig. [l] do not permit a precise determination of the crit- 
ical temperature or any of its associated critical expo- 
nents. The critical point can be more easily found by 
observing the variation of the modulus of the spatial av- 
erage of the field with the temperature: 




where the angular brackets denote ensemble averaging, 
obtained in practice by averaging over many indepen- 
dent measurements, 200 — 1000, corresponding to well 
separated evolution times. The variation of Eq. ([!]) and 
its derivative with temperature is shown in Fig. ^|. In par- 
ticular we see that the derivative diverges at the critical 
point. 
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FIG. 2. The variation of {\<f>\) and its derivative, rescaled by 
1/4, with the temperature. Circles denote measured points. 
The divergence of the derivative marks the critical point. 

The critical temperature can be measured by assuming 
that Eq. (||) behaves like 

(\m\) = A((3-p c ) a , (5) 

which is assumed to hold just below the critical temper- 
ature. By fitting our data at several temperatures to 
Eq. (||), for different choices of the inverse critical tem- 
perature /3 C , we can determine its minimal % 2 value. Si- 
multaneously we obtain the critical exponent a. The 
critical values hence determined are (3 C = 1.906 ± 0.008 
and a — 0.43±0.07. The universal value of a is in reason- 
able agreement with renormalization group calculations 
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Having measured the critical behavior in the fields we 
want to analyze how vortex string densities vary with 
the temperature. We assume ergodicity of the field evo- 
lution once equilibrium is reached. Using this fact we 
analyze at given time intervals the phase of the complex 
field, and associate a vortex to each 2-D lattice cell where 
the U(l) phase winds through 2ir. We then proceed to 
connect the vortices and construct the string network. 
Given a string network, we measure the string density 
in loops and long string, for different values of the bath 
temperature around the phase transition. 

Fig. U shows how the total density per lattice link of 
strings p to t, as well as the density of long and short string 
(pi and pi respectively) changes with temperature. These 
quantities allow us a direct comparison to algorithms of 
cosmic string formation. 




FIG. 3. The lattice link densities of total string ptot, loops 
pi and infinite string pi as a function of inverse temperature 
p . 

We observe a dramatic change in the behavior of the 
string system across the phase transition. The curves of 
fig. H suggest that we can decompose strings into two dis- 
tinct populations, one of loops, say strings smaller than 
a certain cutoff length of the order of the squared 
linear dimension of our total volume and another pop- 
ulation consisting of much longer strings. We will refer 
to the former as the string loop population and to the 
latter, owing to the usual terminology in cosmology, as 
infinite strings. 

At @ c all three densities display abrupt changes in their 
behavior. These changes are phase transitions as in all 
three cases the derivatives of the densities present dis- 
continuities at the critical point. Even more remarkable 
is the fact that the appearance of infinite string seems to 
be linked to the criticality in the fields. The very sharp 
rise in the infinite string density at the critical point may 
suggest that the phase transition associated with it may 
actually be discontinuous in the infinite volume limit. 

In order to clarify this question we measured the infi- 
nite string density variation with j3 for various values of 
Loo. The results are shown in Fig. [I| The temperature at 
which infinite strings first appear clearly depends on the 
choice of L m , getting closer to f3 c as it increases. We also 



verified that for smaller lattice sizes infinite string first 
appeared at lower temperatures, eg. in a 20 3 computa- 
tional domain there was infinite string as soon as = 2.3. 
Both these variations supply evidence that as the volume 
is increased the temperature at which infinite strings first 
appear migrates towards the critical point implying a dis- 
continuous phase transition in pi. At present it is how- 
ever impossible to perform such extrapolation with full 
confidence. 
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FIG. 4. The dependence of pi on inverse Temperature, for 
Loo = 2500, 8000, 20000, 100000 (top to bottom), measured in 
lattice units. 

We can establish that the phase transition in the fields 
and in the string densities coincide by finding the critical 
point and exponents for the latter. 

At the critical temperature the density of strings drops 
suddenly and keeps falling for lower temperatures. The 
decay of ptot which coincides with pi, is excellently de- 
scribed by the exponential law 

PtotiP) = Ptot(Pc)e-^-^ E (6) 

where ptot(Pc) — 0.2 is a universal number (see eg. (^] 
for a very different study yielding the same result) and 
E = 1.836, implying that at temperatures below the crit- 
ical point strings are Boltzmann suppressed. At low tem- 
peratures (in practice for > 6) there is long range order 
and no strings survive. 

Just above the critical temperature the variation of the 
string densities can be best characterized by defining the 
quantity Sp = (p(T) - p{T c )) / p(T c ) . The behavior of the 
infinite string density pi is well described by an ansatz 
similar to Eq. (||), 

5 Pi <x{T-T c yK (7) 

Note the change from f3 in Eq. (||) to T in Eq. (0). We 
find 7, = 0.25±0.01and/3 c = 1.911±0.001, in agreement 
with all previous estimates. For p tot we find: 

5ptotaL{T-T )t**, (8) 

with j3 c = 1.912 ± 0.002 and ltot = 0.39 ± 0.01. 

For temperatures well above the phase transition the 
loop density can be reasonably described by a quadratic 
form 
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p,(/J)~p,(0) 1 
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where pi(Q) — 1/12. This corresponds to a fraction of 
25 ± 1% of all strings in loops for (3 = 0. 

The variations are much more complicated for p tot and 
Pi . The large size of the derivatives close to = make 
it hard to approximate the curves and a successful fit 
can be achieved by a polynomial of very large degree, 
typically 9 or 10, with coefficients of alternating sign. In- 
deed, the behavior looks non-polynomial instead, which 
would mean that there is no standard high-temperature 
expansion for the infinite string density. This is expli- 
cable if we realize that there is no local operator which 
can distinguish an infinite string, and thus we would not 
necessarily expect a polynomial behavior around (3 = 0. 

Particularly interesting is the infinite temperature 
limit. We observe that then the total string density per 
link tends to ptot(0) = 1/3 and that its fraction of long 
strings tends to 75 ± 1%. These results coincide with 
those obtained from the canonical Vachaspati-Vilcnkin 
(VV) algorithms for simulating string formation in the 
early Universe fTo|-|l2|| . Their underlying motivation re- 
sults from the picture that strings should result from ran- 
dom phase fluctuations within spatial domains of a given 
size. Typically, the complex field is given unit modulus, 
with phases assigned randomly to each site on a cubic 
lattice. There are variations on this theme: for example, 
the continuous phase may be approximated by three val- 
ues jl0| , which leads to a slightly different string density 
and infinite string fraction. Nonetheless, assigning phases 
randomly to each lattice site corresponds precisely to the 
infinite temperature limit, and thus VV algorithms gen- 
erate ensembles of string at infinite temperature. These 
ensembles are used as initial conditions for the free evolu- 
tion of the string network, using a numerical solution to 
the classical dynamics of relativistic strings [[!) . Thus nu- 
merical simulations of this type model an instantaneous 
quench from infinite to zero temperature. While these 
initial conditions are ultimately unphysical, the network 
of strings is observed to approach rapidly a self-similar 
evolution, which serves to disguise the initial state. On 
the other hand, if one would like to describe a network of 
strings below the critical point the Vachaspati-Vilcnkin 
algorithm would be a completely inappropriate starting 
point. 

Finally, our observation that the phase transition in 
the fields may coincide with the point when infinite string 
starts forming, suggests that strings do indeed play a 
fundamental role in the critical behavior of the system, 
as has been suggested so often in the literature [[13] [L6| . 
On the other hand, statistical studies of string networks 
Ji~7| commonly identify the critical point in string den- 
sities with the Ginzburg temperature Tq < T c . Above 
we presented evidence for the coincidence of the critical 
behavior in the fields and in the strings at T c . In a fu- 
ture publication [j^J, a closer look at the properties of 
the string length distribution and a more detailed scale 



dependence study will help to clarify these problems. 
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